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We present a new method for locating unstable periodic points of one dimensional chaotic 
maps. This method is based on order statistics. The densities of various maxima of the 
iterates are discontinuous exactly at unstable periodic points of the map. This is illustrated 
using logistic map where densities corresponding to a small number of iterates have been 
obtained in closed form. This technique can be applied to the class of continuous time 
systems where the successive maxima of the time series behave as if they were generated 
from a unimodal map. This is demonstrated using Lorenz model. 

PACS numbers: 05.45. +b, 02.50.-r 
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Existence of a dense set of unstable periodic orbits 
(UPO) is one of the characteristic properties of a chaotic 
system jjj. These orbits represent the skeleton for the 
strange attractor of dissipative dynamical systems. Many 
quantities that characterize chaos in the system, such as 
the fractal dimension, average Lyapunov exponent, en- 
tropy and the invariant measure of the corresponding at- 
tractor can be determined by knowing the properties of 
UPO @§J|||§]. Extraction of UPO is a necessary step 
in several studies. For example, knowledge of the loca- 
tions of various cycles is necessary for control of chaos 

0. Cycles are found to be useful in the synchronization 
of chaotic signals ||. Moreover, the quantum mechani- 
cal properties of classically chaotic conservative systems 
have, in the semiclassical regime, a series expansion with 
respect to the lengths and stability coefficients of the pe- 
riodic orbits The importance of unstable periodic 
orbits and their detection has attracted the attention of 
several researchers and a number of numerical methods 
have since been developed to extract unstable periodic 
orbits g0,|||||l|Jl|Il. 

In this letter we present a new technique to extract 
the UPO of one dimensional chaotic maps using order 
statistics. This method works for all one dimensional 
chaotic maps for which the existence of invariant density 
is guaranteed. 

The theory of extreme values and its generalization to 
order statistics is a classic subject and is extensively used 
in the study of independent and identically distributed 
random variables |Q. Let {Xi, Xs, - ■ • , X n } be an n- 
point set and M% (1 < k < n) be its k th maximum 

1. e., M™ < Af™- 1 -.. <Ml<---Ml< M*. Order 
statistics is the study of distributional properties of 

In the context one dimensional chaotic map, it has been 
shown that the extreme value density is discontinuous on 
a set of points belonging to its unstable periodic orbits 
p7f. However, this method cannot locate all the periodic 



points. This has motivated us to investigate whether all 
the unstable periodic points can be extracted by using 
order statistics instead of simple extreme value analysis. 

Let f(x) : [a,b] — > [a, b] be a continuous, one dimen- 
sional chaotic map with an invariant density p s (x). Let 
{x Ql xi = f (x ) , • • • ,x n -i = (x Q )} be an n-point 
set. The density of the k th maximum is denoted by p„(x). 

Let be the set of locations of the discontinuities of 
p k n . We observe that O n = \J% =1 S* for all n, where O n is 
the set of interior (excluding the end points a,b) periodic 
orbits of all orders strictly less than n. An outline of the 
proof of the above result is given towards the end, see 
for details . This remarkable correspondence between 
the densities of the order statistics and UPO enables one 
to extract the periodic points of any order by choosing 
an appropriate n. The locations of the cycles are, thus, 
read off from the extremely sharp discontinuities of the 
densities. 

We now illustrate the method by applying it to logistic 
map, x n+ i = f(x n ) : [0,1] -> [0,1] = \x n (l - x n ), for 
A = 4, a well studied unimodal map exhibiting chaos [0. 

We analytically derive the expression for the densities 
of the first and second maxima of a three point set of 
logistic map using the formulation described in equation 
(|). They are 
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Note that the point x = (5 — v5)/8 is not present in 
S3 = {3/4, (5 + \/5)/8}, but is a periodic point of order 
two, which clearly shows that extreme value statistics is 
not adequate to extract all the UPOs and one needs order 
statistics. 

Numerically, given a map with an invariant density (a 
mere existence is necessary and sufficient to render mean- 
ing to the averaging procedure), one computes in the 
following way. Starting with an initial condition, n — 1 
successive iterates of the map are obtained and arranged 
in the ascending order of their magnitude. The k th max- 
imum is thus picked up. This is repeated for many initial 
conditions to obtain a histogram representing p\. For 



the logistic map p 
and is shown in fig 
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FIG. 1. Numerically computed p\ for k = 1,2,3 and 4 
of the logistic map (10 7 initial conditions, bin width=0.001). 
The locations of the discontinuities and the cycle points of all 
orders < 4 are compared. The cycle points are indicated by 
vertical lines. Note the absence of discontinuity at the origin 
which is a fixed point (see the text). 

Further, the topological entropy, which is a quantita- 
tive characterization of chaos in the system, can be esti- 
mated |]13] from the number of periodic orbits as follows, 
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where N(n) is the total number of periodic points of or- 
der n. We denote the number of elements in a set A 
by Card(A). Since N(n) = Card(O n+ i) - Card(O n ), it 
can be obtained from the discontinuities of p t ^ l+1 {x){k = 
1 • • • n + 1) and that of p^(x)(k = 1 • • • n). In the case of 
logistic map we obtain Card(O n ) as 3, 9, 21, 51 and 105 
for n =3, 4, 5, 6 and 7 respectively. Since this formalism 
is based on the discontinuities of the density over an in- 
terval [a, 6], it cannot indicate a periodic point occuring 
at the extreme points of the interval a and b. Thus, in 



the above list (see also fig [j]), the fixed point at the ori- 
gin is not included. The topological entropy calculated 
is 0.6648, which is within 4% of the exact value ln2. 

However, for A < 4, the invariant density, p s (x) is not 
available in closed form. The existence of p s {x) can be 
checked and is sufficient for computing order densities. 
For most A < 4, p s (x) itself is discontinuous on a large 
number of points. Thus, the corresponding order densi- 
ties pick up discontinuities at those points at which p s {x) 
is discontinuous and also at unstable periodic points. 
Since the strength of the inherent discontinuities of p s (x) 
in Pn(x) is large compared to that of the discontinuities 
at periodic points, it is numerically difficult to count or 
detect their locations. 

There are situations where the essence of the dynamics 
of a continuous time dynamical system is captured effec- 
tively by one dimensional maps or their equivalent. In 
his classic paper Lorenz showed that the successive peaks 
of a one dimensional time series behave like iterates of a 
map pel , see also [H . We work with the Lorenz model 



x = a (y — x) , y = x (r — z) — y, z = xy — bz (4) 

for parameters r = 28, a = 10 and b — 8/3. 

We compute the p\ (x) for different n for the map con- 
structed out of the successive maxima of the time series 
corresponding to the state space variable z. Existence of 
an invariant density for the so constructed map is veri- 
fied numerically. p\(x) of this map also is discontinuous 
on a set of points. The topological entropy of the map is 
calculated using equation (|J). For n=3, 4, 5 and 6, the 
topological entropy /i=0.6486, 0.7361, 0.7700 and 0.7708 
respectively, showing a reasonable convergence in the nu- 
merical value of topological entropy. 

It is important to consider the effect of noise on the or- 
der densities to explore the possible applicability of this 
method to an experimental time series convoluted with 
noise. Here, the distribution and strength of the noise 
play a decisive role. Our preliminary investigations in- 
volve addition of noise, generated from uniform distribu- 
tion, to the iterates of the logistic map with A — 4. The 
strength of the noise is 4% of the variance of the iter- 
ates. The spill over from [0,1] is reinjected. The strength 
of the discontinuity at the periodic points is observed to 
diminish with increase in the strength of the noise. 

We now turn to establish the contribution of periodic 
points to the order statistics rigorously. The formulation 
of the order statistics involves the probability P% (x) that 
e (x, x + dx) and is given by 



(k 



1)! ^ ^ 



X < Xi t < X 



Prob 



dx 

i > x\i el 
< x\j e J 



(5) 



2 



where X = {t2,---,ik} and J = {i k+1 , ■ ■ ■ ,i n } with 
{ii, • • • , i n } being a permutation of {0, • • • ,n — 1}. Us- 
ing the joint probability density of the n points together 
with the above equation, we derive an expression for the 
cumulative distribution F k (x) of the k th maximum . The 
order density p k (x) is the derivative of F k (x). 

Our central result is : given a continuous one dimen- 
sional map / : [a, b] — > [a, b] with a continuous invariant 
density p s (x), the order density p k (x) is discontinuous 
at an interior point if and only if ir* is an unstable 
periodic point such that f l {x±) — x* with I < n — 1. 

We give an outline of the proof and the notation used is 
described here, (i) A periodic point of order I is denoted 
by pjip, i.e., fipjip) = p 3 ip. The index (3 = 1 ■ • ■ I corre- 
sponds to the distinct points of the Z-cycle. In a chaotic 
map, in general, there exist more than one Z-cycles, and 
the number of /-cycles JV; increases with I. The index j 
runs over the different Z-cycles and hence j — 1, • • • , JVj. 
For example, in logistic map, we have two orbits of pe- 
riod 3. The periodic points of a given /-cycle are ordered 
such that pju = maxjpj/^} corresponds to the maximum 

of that cycle, (ii) Let y = p l {x). The set of preimages 
of x with respect to p are {gi a (x)\P {gia (x)) = x\ ,a = 
l,---,7j where Ii is the number of preimages. For ex- 
ample, the preimages of x = 3/4 of the logistic map are 
{3/4,1/4}. 

As a first step, we show that the extreme value den- 
sity, Pn(x) is discontinuous at x if and only if x is 
the maximum point of a periodic orbit, i.e, x = Pju, 
I < n — l,j = 1, • • • JVj. Proving this involves two stages 
namely, to show that the extreme value density cannot 
be discontinuous at x unless x is the maximum of a pe- 
riodic orbit and then show that is discontinuous 
at the maximum of all the periodic points. The extreme 
value density can be formulated by using equation (|^) to 
be 

n— 1 n— 1 

Pl{x) = Ps(x) Y[e{x-f (X)) + E E A ^ x ) 

j—1 i—1 a—1 

9 (x - g ia (x)) l[e(x- P (g ia (:r))) (6) 

where A ia (x) = p s (g ia {x)) /\-^p{x)\ x=g ^ [x) . The 
above equation can be recast as 

pl(x) = Ps {x)C {x) + E A ia {x)C ia {x) (7) 

i—1 a — 1 

It can be shown that if p s (x) is continuous then 
Ai a (x) is also continuous, see |lj| for details. As all 
the terms in equation (^|) are positive definite, proving 
either Cq (x) or Ci a (x) is discontinuous at a point is suf- 
ficient to prove Pn(x) to be discontinuous. Also note 
0(z) = (1 + sgn{z))/2 is discontinuous only at z = 0. 

Since p s (x) and hence Ai a (x) are continuous, Pn(x) is 
discontinuous only if either (a; — y (x)) for some / < 



n— 1, or (a; — gi a (x)), or (x — f m (gi a (x))) for some 
m ^ i,m < n — 1 is discontinuous . 

If (x — p (x)) is discontinuous then x = f l (x) im- 
plying x S {Pjip}- Considering Co(pjip + e) (where e is 
an infinitesimal quantity of appropriate sign), we have 



Cofav + e) = [] e (PUP ~ f" (PW) + 0(e)) (8) 
fe=i 

If ^7/3 7^ Pj-jj, in the product, there exists a fej such that 
f kl {pjip) = Pju and the product vanishes as pjj; > pjip 
for all /? < I. This shows that Cq(x) cannot be discon- 
tinuous unless x = pju . 

Q(x — gi a (x)) is discontinuous at x = gi a (x). We 
have p(x) = p (gi a (x)) = x and this implies x G {Pjip}- 
Considering Ci a (pjip + e), the leading order term is 



C ia (pji{3 + e) = 



1 dx® ia I k= Pj*£ 



n e (^- /* fotf) + coo) w 

It can be shown that 3i Q (a;) = 1/ (#)) and 
|/'(5ia(Pji/3))| is always greater than one for all the un- 
stable periodic orbits. This means the first function 
in the above equation is always non zero. Thus, by a 
similar reasoning outlined above, the second term in the 
product cannot survive unless x — pju . 

Consider (x — f m (gi a (x))) which is discontinuous 
at x = P n (gi a (x)) for some m e {1, ■ • • ,n — 1}\{«}- 
There exists an I such that x = f m (gi a (x)) implies x = 
f (x), where / = max{m — i, i — m} implying x £ {Pjip}- 
By using the previous arguments, one can show that the 
discontinuity cannot occur unless x = pju. This proves 
that p\(x) cannot be discontinuous unless x — pju for 
some i&j. 

Conversely, we show that at every x = Pjiu Pn{x) is 
discontinuous. This is done in two steps for convenience, 
namely (a) if pju is such that I > (n — l)/2, Cq(x) regis- 
ters a discontinuity at x — pju while (b) if I < (n — l)/2, 
Ci a (x) will be discontinuous at x = pju . 

(a) Consider Cq (pju + e) and writing it as a product of 
three terms, we have 



Cq(p 3 u 
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\{®{p 3 ii- f k {Pju) + 0(e)) (10) 

n-l-1 

II ®{v 3 u-f k+l {p ]U ) + 0{e)) 



fe=i 



f k (Pjll ) < Pjll f° r all fc < Z — 1 and n — I — 1 < I — 1 if 
I > (n— l)/2. This implies that the product terms in the 
equation dTfl) are non zero and the first term also can be 
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made non zero by choosing an e of appropriate sign. This 
proves that if I > (n — l)/2, at x = Pjzi,^ < n — 1, j = 
1, • • • , Ni, Cq(x) and hence Pn(x) are discontinuous. 

(b) We prove that, given x = Pj mm ,'m 7^ i and m < 
(ri— 1)/2, there exists a function in d a (x), see equation 
(||), with argument x — f l (gi a (x)) , I ^ i such that I — 
i ± m which is discontinuos at 2 = Pj mm , see for details 
p8[ implying Ci Q (a;) and hence p^(a^) to be discontinuous 

at X — Pjmm- 

Similarly one can prove that the density corresponding 
to k = n, the minimum of the iterates, Pn{x) is discon- 
tinuous at x if and only if x — Pjn, i.e, at the minimum 
of the periodic orbits. 

Finally, it can be proven that at every periodic point 
x = Pjifj, there exists a k such that the order density 
p„{x) is discontinuous at Pjip- This result is obtained by 
subjecting p^ix) to a similar analysis outlined in the case 
of Pn(x), see Q for details. 

In summary, we have presented a new method to ex- 
tract all the unstable periodic points of one dimensional 
chaotic maps. This method uses order statistics. We cal- 
culate p\ of the logistic map analytically for n = 3 and 
numerically for n > 3 and illustrate that their disconti- 
nuities coincide with the unstable periodic points of all 
orders less than n. We demonstrate the applicability of 
this method to the map constructed out of the successive 
maxima of the time series z(t) of Lorenz equations. In 
both of the above cases topological entropy has been es- 
timated. Further, we give an outline of the proof which 
establishes this connection between order statistics and 
unstable periodic orbits for any continuous map with a 
smooth invariant density. 
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